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Abstract 



After a brief review of the present status of nonextensive statistical me- 
chanics, we present a conjectural scenario where mixing (characterized by the 
entropic index q m ix < 1) and equilibration (characterized by the entropic in- 
dex q eq > 1) play central and inter-related roles, and appear to determine a 
priori the values of the relevant indices of the formalism. Boltzmann-Gibbs 
statistical mechanics is recovered as the q m %x = Qeq = 1 particular case. 

Human knowledge progresses along very many paths, and very rarely these paths follow 
a systematic and logical ordering. What we can presently witness about the formalism 
frequently referred to as nonextensive statistical mechanics (first proposed in 1988 ||T| and 
further implemented in |2],|3|; for reviews, see |§-[§]) is by no means exception. Since it is 
the purpose of the present lines to start with a review, let us present it in an order which 
makes some epistemological sense, although it does not necessarily follow the chronology 
of the events. We will successively comment on (i) mathematical properties, which make 
the formalism to appear just as applied mathematics; (ii) the successful confrontations with 
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experimental (as well as computational) results, which we believe raises the formalism to the 
status of theoretical physics, since it has a compromise with phenomena indeed occurring in 
nature; and (hi) the lines along which the calculation of the relevant entropic indices q m %x 
and q eq appears as possible, so that the formalism becomes a closed theory (in other words, 
a complete theory: see [[?]]), where everything can in principle be calculated a priori once the 
dynamics of the system is fully characterized. 
Applied mathematics 

The axiomatic starting point is the proposal of a possible generalization of Boltzmann- 
Gibbs (BG) statistical mechanics by postulating the following entropic form 

1 _ yw n i i w 

s q = fl y =H-> (E» = i). (!) 

q i Pi i= x 

where we are using ks = 1 (without loss of generality), q G 1Z, (■ ■ ■) = Y^LiPii' ") an d 
ln ? x = xl - L ^~ 1 (with hi! a; = lnx) [consistently, its inverse function is = [1 + (1 — q)x] 1 ^ 1 ~ q \ 
with = e x \. Clearly, for q — 1, we recover the usual expression Si = — X^iPilnpi, from 
now on referred to, for simplicity, as the BG entropy. Also, if A and B are two probabilis- 
tically independent systems (i.e., pfj +B = pfpf), S q satisfies the following pseudoadditivity 
property: S q (A+B) = S q (A) + S q (B) + (l-q)S q (A)S q (B), hence q = 1, > 1, < 1 respectively 
correspond to extensive, subextensive and superextensive cases. 

For an isolated system (microcanonical ensemble), optimization of S q yields equiproba- 
bility, i.e., Pi = 1/W (Vz), hence 

S q = ln q W . (2) 

Before going on, it is important to stress that the pseudoadditivity property mentioned 
above does not exclude that, for some classes of interdependency between A and B, a special 
value q* could exist such that extensivity could be recovered in the sense that S q * (A + B) = 
S q *(A) + S q *(B) . If so, we could say that the adequacy of the entropic form enables the 
preservation of the traditional property of extensivity for the entropy. As mathematical 
illustrations of this fact, let us consider two possible cases at equiprobability. First, if 
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W ~ afi N (with a > 0, \i > 1 and iV = number of elements — > oo), then Si = lnVT ~ 
(ln//)iV oc N, hence g* = 1. Second, if W ~ 6iV (with 6 > 0, z/ > and iV -> oo), 
then, for g < 1, S q = W ^ q ~ l ~ ~ ^A^ 1 "^, consequently the choice q* = 1 - l/v 
implies S q * oc N, once again recovering the traditional proportionality between entropy and 
N, when iV >> 1. 

For a system in thermal equilibrium with a thermostat (canonical ensemble), we must 
optimize S q with a further restriction, namely |§ 



J2 p ^ = u q , (3) 

i=l 

where the escort distribution Pi = p\j Y^=iP% { e i} are the eigenvalues of the corresponding 
Hamiltonian (with the associated boundary conditions), and U q is a fixed finite value for the 
generalized internal energy. The optimizing distribution is then given by 

e -(ei-U q )/T q 

Pi = sriv -te-u q )/T q ' ( 4 ) 

where T q = TY^LiP q j , l/T = [3 being the Lagrange parameter associated with constraint 
(3). We can verify that, for all values of q, Pi is invariant under shifts of the zero of the 
spectrum of energies {ej}. Also, for q — 1, we recover the usual BG equilibrium distribution 
Pi = e~ e ^ T I X)J=i e~ tj / T . Eq. (4) can be rewritten in another form, namely 

e~ q llT ' q 

Pi = J _ e . /T / > ( 5 ) 

with T' q = T g + (1 — q)U q . For fixed and large values of T' q , we have pi ~ (1 — ei/T' q )/\Y^ =x {l — 
6j/T')] for all values of q. In other words, all nonextensive equilibrium thermostatistics share, 
at large temperatures, a common distribution which can be equivalently considered to be 
the BG one §. 

These are the essential steps. From these, many other have been developed for arbitrary 
q, such as the if -theorem, Ehrenfest theorem, Bogolyubov inequality, factorization of the 
likelihood function, Onsager's reciprocity theorem, Kramers and Wannier relation, Pesin 
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theorem (conjecture), fluctuation-dissipation theorem, and others. Also, a variety of tra- 
ditional statistical-mechanical techniques for treating many-body systems is now available 
for arbitrary q, such as Green functions, variational method, perturbative methods, path 
integrals, Lie- Trotter formula, simulated annealing, and others. For review of all these, see 
and references therein. 
Finally, in Table 1, we present the most relevant historical steps of the foundations 
of BG statistical mechanics, and their generalizations for arbitrary q. In 1860, Maxwell 



presented its celebrated Gaussian distribution of velocities ||1 1|| . In 1872, Boltzmann [12 
arrived, through the molecular chaos hypothesis (stosszahlansatz) , to the celebrated expo- 
nential weight as the stationary state of his partial derivative kinetic equation for distribu- 
tions in the so called /i-space (projection, on the one-particle phase space, of the states of all 
particles). In 1902, Gibbs [H| presented how, within a variational principle using entropy 



as the relevant functional, the exponential weight can be reobtained, this time in a more 
general framework, namely in the so called T-space (phase space of all particles). Gibbs 
equilibrium distribution was later reobtained through a variety of other manners, namely 



by Darwin and Fowler in 1922 flTJJ using a steepest descent argument, by Khinchin in 1949 



|T5| using the law of large numbers, by Balian and Balazs in 1987 |T(| and by Kubo et 



al in 1988 fT7fl, performing countings in the microcanonical ensemble (isolated system). In 



parallel with these developments, Shannon in 1948 |18| and Khinchin in 1953 estab- 
lished necessary and sufficient conditions for the entropic functional to be — X^iPjhi p^ 
All these arguments have been generalized for arbitrary q. Gibbs path was followed in 1988 
and thereafter |l]-|3|], the Darwin- Fowler, Khinchin and Balian-Balazs paths were followed 



in 2000 by Abe and Rajagopal in [20], |^TJ and |22| respectively, and the Kubo path was 



followed by Abe and Rajagopal in 2001 |23| . Boltzmann path was followed in 2001 by Lima, 
Silva and Plastino |M] as well as by Kaniadakis in ||25|| . Shannon and Khinchin paths for 
the necessary and sufficient for the entropic form were respectively followed by Santos in 



1997 [25| and by Abe in 2000 [27 1. All these generalizations, without exception, proved to 



be consistent among them and consistent with Eqs. (1-4). 



At this point it seems appropriate to comment that, as a whole, the formalism exhibits 
a remarkable mathematical "texture" , which quite naturally extends to arbitrary values of 
q the properties that are since long known for q — 1. 
Theoretical physics 

Let us now argue in the sense of transforming the above mathematical formalism into 
theoretical physics by comparing theoretical with experimental results. We must however 
warn the reader that many of the available applications of nonextensive statistical mechanics 
concern open systems (stationary states of open systems), and not only (meta) equilibrium 
states of the time-independent Hamiltonian systems addressed in Eqs. (3-5). To be more 
precise, all thermodynamical equilibrium states are stationary solutions of some family of 
partial derivative equations (e.g., of the Boltzmann kinetic equation, for classical systems), 
but the opposite is not true. There are stationary states which are not tractable in thermo- 
statistical terms, this is to say susceptible of being founded in geometrical aspects of some 
phase space (or of some Riemann subspace in that phase space) or analogous spaces such 
as the Hilbert or Fock ones. To illustrate what we mean about geometrical founding, it 
is certainly instructive to recall that it is the symmetrization or the antisymmetrization of 
the corresponding wave functions that determines the transmutation of M axwell-Boltzmann 
statistics into Bose- Einstein and Fermi-Dirac ones, consistently leading to entropic function- 
al which differ from the classical one. The point is that if the dynamics of the system is 
completely known, its possible stationary states always are in principle calculable, but for 
some classes of such systems it is not necessary to follow its dynamics: we can directly, 
geometrically, calculate the stationary state, which can then be considered as optimizing 
some entropic functional. Statistical mechanics focuses on such states. The difficulty is 
of course to know a priori what specific statistics is to be applied to a given system. We 
come back onto this point later on. But at the present stage, let us emphasize that in our 
understanding statistical mechanics emerges if and only if dynamics can in some sense be 
replaced by geometry. Unfortunately it is by no means clear when this is possible, but we 
shall refer to it as geometrizable dynamics. As an attempt to clarify these arguments, we 
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schematically display in Figs. 1 and 2 respectively the traditional and the present views 
on the connections between dynamics and statistics (see also Cohen's contribution to the 
present proceedings). In Fig. 1 it is stressed that, in the traditional view, an unique type 
of thermal equilibrium exists and this is the BG one (either in the micro-canonical, the 
canonical or the various grand-canonical forms). We believe that statistical mechanics is 
wider than that, as indicated in Fig. 2. Every time that dynamics (of a finite or infinite 
system) can be automatically taken into account by geometrical considerations, theory of 
probabilities can, for a variety of purposes, efficiently replace the knowledge of the time 
evolution of the mechanical system. Within classical BG statistical mechanics, dynamics of 
vast classes of isolated systems (typically involving short-range interactions) can be replaced 
by the hypothesis of equiprobability in the occupation of the accessible phase space together 
with a connection of macroscopic entropy with the relevant phase space volume. We are of 
course referring to the Boltzmann principle S — k In W (see also Gross' contribution to the 
present proceedings). It is however physically appealing to think that more subtle situations 
can also be handled within statistical-mechanical methods. Such could be the case when, in 
addition to the knowledge of the volume W, we need to characterize a physical, dynamical 
bias : this is the role of q within nonextensive statistical mechanics. The just mentioned 
Boltzmann principle would have to be generalized in such cases by S q = k\W l ~ q — 1]/[1 — q\- 
Of course, there is no reason for thinking that such geometrization of dynamics could not 
be in principle done for other, possibly more complex, systems, outside of the g-statistics 
focused here. Within this scenario, nonextensive statistical mechanics appears to be just the 
first non-Boltzmannian thermostatistical formalism; many others are in principle thinkable, 
corresponding to various manners for replacing (specific classes of anomalous) dynamics by 
geometry. In our present formalism it is clear that W is a geometrical concept, but the 
reader might be puzzled by the fact that we are including q in the same category. This point 
will become transparent soon, when we shall illustrate how multifractal geometrical consid- 
erations enable us, at least for some simple specific cases, to uniquely determine q a priori 
from the mechanical characterization of the system (and not only from fitting procedures, 
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as frequently done, faute de mieux, in the experimental literature). 

To make connection with nature, it is mandatory to mention that the present formalism 
has been successfully applied to a considerable variety of systems, such as Levy [^] and 
correlated J^J anomalous diffusions, peculiar velocities in spiral galaxies |ftm , turbulence in 
electron plasma [[£[[], fully developed turbulence |3l^P5], citations of scientific papers [|3f) 



linguistics |37 |, reassociation in folded proteins f38 |, quantum entanglement [|3^,f|(|, electron- 
positron annihilation [41]], quark-gluon plasma [ 42| , cosmology [^3j^4| , hadronic scattering 
f45fl , motion of Hydra viridissima ||46|| , low- dimensional maps []47l , f48|l , inertial classical planar 



rotators ferromagnetically coupled at long distances |]49| , p0|| , among others. To be more pre- 
cise, this formalism addresses systems which, in one way or another, exhibit some relevant 
multifractal structure. This can occur through a variety of physical mechanisms, such as spa- 
tial and/or temporal long-range interactions, mesoscopic dissipation, multifractal boundary 
conditions, quantum entanglement, and others. 

This is a good point for warning the reader about the fact that stretched exponentials 
and q- exponentials can be numerically very close to each other for intermediate values of 
the abscissa (see Fig. 3), although they are definitively very different both close to the 
origin and approaching infinity. It is important to realize these features in order to really 
appreciate the strength of experimental and computational evidences favoring one or the 
other functional form whenever fittings are involved. The present formalism naturally leads 
to g-exponentials, rather than to stretched exponentials (frequently used for fittings during 
the last 10-15 years). But only fittings that are satisfactory over relatively large physical 
ranges can be acceptable in order to distinguish between these two functional forms. For 
example, in the case of Hydra viridissima just mentioned | 46fl , both a g-exponential and a 
stretched exponential fit well the experimental distribution of velocities if only relatively 
small velocities are taken into account. In this particular case, it is because large velocities 
were experimentally measured as well, that it became possible to dismiss the stretched 
exponential function, and retain the g-exponential one. 
Closed theory 



We shall from now on note q eq the value of q characterizing the equilibrium distribution 
(or the stationary distribution if the system is open); typically q eq > 1 (see, for instance, 
p^ , |35|^T| - p|p6| . pO[ ) . Let us also address a different quantity, noted q m ix, related to the 



mixing properties of the system; typically q m i X < 1 (see, for instance, p3| , ^7| ). The relation- 
ship between these two different values of q is under intensive study nowadays (see |]51|| ). 
Our focus here is, as mentioned previously, to develop a scenario where g m j X , q eq and the 
metaequilibrium states that frequently emerge for some nonextensive systems, play deeply 
entangled roles. The final outcome is to illustrate the lines along which dynamics can a 
priori determine the values of q to be used for specific physical models, so that the whole 
thermostatistical formalism becomes a closed and complete theory. 

To start uncovering the scenario it is enough to consider systems whose phase space 
(noted x) is one-dimensional, e.g., one- dimensional maps such as the logistic one. If we 
note Ax(t) the discrepancy, as a function of time, of two trajectories initially discrepant 
of Ax(0), we can define the sensitivity to the initial conditions (or mixing function) £(£) = 
liniAa:(o)-<-o Ax(t)/ Ax(0). The most frequent case is that where £ satisfies the following 
differential equation: 

where Ai is the Lyapunov exponent (the subscript 1 will become clear soon). It follows that 

e = e Al ' (7) 

Positive and negative values for Ai respectively correspond to strong sensitivity and insensi- 
tivity to the initial conditions. What happens when Ai = 0? The generic answer is that Eq. 
(6) is not applicable anymore, and we must take into account the next term, so we consider 
now 

f = K*. Z qmix > ( 8 ) 

where we focus on the case \ qmix > and q mix < 1. The solution now is 



e = [i + (i - ft^) ^ mi ,t\ 1/{1 - qmi ^ , (9) 

which reproduces solution (7) if q m i X = 1, and asymptotically behaves like tVl 1 -?™.™) if 
Qmix < 1- The latter will be referred to as weak sensitivity to the initial conditions (or weak 
mixing). Let us finally consider the most general case along this line, namely 

§ = Ai£ + (A f ,_-A 1 )e"% (10) 



whose solution is 



£ _ ^ Qmix _|_ Itnix g(l—q m ix))>lt 1 ~1mix (H) 

L Ai Ai i 
This function is illustrated in Fig. 4 in such a way as to exhibit the crossover from the 
power-law regime at intermediate times to the exponential regime at long times, which 
occurs when < Ai << X qmix and q mix < 1. More precisely, for times t satisfying < 
t « t* = 1/[(1 - q m ix) \ mix ], we have an integrable-like regime (with £ ~ 1 + X qmix t) 
(i.e, characterized by q = 0), for times satisfying t* « t « t** = 1/[(1 — q m i X ) X\] 
we have a power-law regime (with £ ~ [(1 — q m i X ) X qmix t] 1 ^ 1 ^ qmix \ i.e., characterized by 
Q — Qmix), and finally for times satisfying t » t** we have an exponential regime (with 
£ x irnix e Ai^ j e ^ characterized by q = 1). These facts lead to the following nonuniform 
convergence: linr^-o lim^oo £f = lim Al ^ lim^o = oo, whereas lim^oo lim^-o 



lim^oo —1 — ^4 = —1 — < oo. Analogously, using the g-logarithm function ln„x, we have 

-L Qmix t 1 Qmix 

lim Al ^ lim^oo ln<? "p ^ = lim Al ^ lim^oo ^'/(i-g™*) = oo, whereas lim^^ lim Al ^ ^f 2 ^ = 
^qmix < oo. In other words, for 1 << t « £**, we have - Sx f s - ~ X qmix , whereas, for t » t**, 
we have ^ ~ Ai. We believe that these features constitute the basic scenario of validity of 
Boltzmann-Gibbs statistical mechanics versus validity of nonextensive statistical mechanics. 
This belief is supported by at least three examples, namely the standard map (see and 



references therein), another, billiard- inspired, two-dimensional conservative map and 
the system of N classical inertial planar rotators coupled all with all through long-range 
interactions (so called a-XY model, which for a = reproduces the HMF model; see P9]J5H 
and references therein). 



The standard map is a conservative map with a two-dimensional phase space (the lowest 
dimension at which a map can be conservative). It includes a nonlinear term introduced by 
a constant a. In the limit a = the system is integrable, and for a 7^ it is chaotic (i.e., it 
has over entire regions of the phase space two nonzero Lyapunov exponents of opposite sign 
and equal in absolute value). The positive Lyapunov exponent Ai is a monotonic function 



of a which vanishes for a = 0. In the limit < a « 1, this map precisely exhibits |52 
the scenario described above with q mix ~ 0.3 as studied through the time evolution of the 
entropic form S q . Indeed, for times below a crossover time (which appears to diverge when 
a approaches zero), 5*0.3 increases linearly with time (see |52| for details), whereas for times 
above that crossover time, it is S\ which linearly increases with time. The same behavior is 
observed for the above mentioned billiard-like map, but with q mix ~ 0.5. 

Let us now turn onto our third example. The system of rotators mentioned above has 
been studied numerically and presents, in the microcanonical ensemble (i.e., isolated) a 
second order phase transition at some total energy (conveniently scaled with N). This 
system presents anomalies both above and below that critical energy. Above that point, the 
entire Lyapunov spectrum collapses to zero when N — > 00. This is to say, 1/N plays a role 
analogous to a in the map just described. Below that critical point, the system exhibits at 
least two (probably only two) basins of attraction with respect to the initial conditions at 
which the system is dynamically started. There is a basin of attraction which exhibits an 
equilibrium at the temperature that the recipe of BG statistical mechanics provides, and 
whose distribution of velocities precisely is the Maxwellian one. But if we start from the other 
basin of attraction, the system equilibrates at a finite temperature different (below) from 
that indicated by BG statistics, and its distribution of velocities exhibit unusual long tails. 
After some long time the system crosses over essentially to the BG equilibrium state. The 
duration of this anomalous metastable state diverges with N, in a manner which, once again, 
strongly reminds the scenario we advanced in the present paper. Indeed, if we consider the 
limjv^oo lim^oo case, the equilibrium appears to be correctly described within BG statistical 
mechanics, but if we consider the lim^oo liniAr^oo case, this is definitively not true, and 
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some other thermostatistical description becomes necessary, apparently the nonextensive 
one. To make the scenario stronger, it can be verified that, during this metastable state, 
the Lyapunov spectrum also vanishes in the limit N — > oo. 

We have sketched above what happens with mixing. What can we say about equilibration 
and q eq l In other words, how physical quantities relax onto the corresponding equilibrium 
values? Let us refer once again to the logistic map as an illustration. In the region where 
strong chaos exists (hence above the chaos threshold), for fixed W and using an ensemble 
where all initial conditions belong to a single among the W windows of the partition, S\ (t) 
exhibits a linear increase with time and then a saturation at S\(oo). Moreover we verify 
that o\ = |1 — Si(^>) I ~ e~*/ T1 , T\ being of the order of 1/Ai (consistently with Krylov's ideas 
more than half a century ago [54|). In other words, 0\ essentially satisfies ^ = —<J\/t\. 



s (t) 

However, at the chaos threshold, Ti diverges, and the time evolution of a a . = 1 — g qmi , \ 

' J J. O J Iran I 6 - (OO) I 

is essentially given |>lj by d<Tq ™ ix = -(v Qmix ) qeq /r qeq , with q eq > 1 and r qeq hopefully of the 



order of 1/X qmix . The solution of this differential equation is given by a qmix = 1/[1 + (q eq — 
1) ~k/ T qeq\ i which reproduces e~*/ Tl for q eq = 1. Furthermore, q eq depends on W, and, 

in the W — > oo limit (infinitely fine graining), we observe (within some degree of accuracy) 



5I[ the following finite size scaling 



g eg (oo) - q eq {W) oc ^- , (12) 

with q eq (oo) > 1. This fascinating relation has up to now been verified only for the ^-logistic 
maps. For these maps and all values of z that have been checked, the values for q eq (oo) 
precisely coincide - supreme suggestion of correctness of the present conjectural scenario! 



with the values numerically obtained in [48|, where the initial conditions were spread 



uniformly all over the entire accessible phase space (in a typically Gibbsian manner). 

We are unfortunately not in position to rigorously prove the validity of the above rela- 
tions nor discuss the detailed hypothesis they must involve. It has however been possible to 
give some physical consistency to the conjectural scenario by analyzing three different phe- 
nomena, namely related to electron-positron annihilation [[y]], to fully developed turbulence 
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35H and to the Henon-Heiles Hamiltonian [04 . 



The distributions of transverse momenta of hadronic jets produced by electron-positron 
annihilation experiments at CERN have been discussed by Bediaga et al in |^T[. The values 
for q eq they obtained depend on the collision center-of-mass energy E, which plays a role 
similar to W since the graining is finer for increasingly high energies. The data in their 
Table 1 can be organized as indicated in the present Fig. 5, with q eq (oo) ~ 1.30, and the 
role played by the exponent 1/2 is that of q m i X in Eq. (12). 

The distributions of velocity differences in fully developed turbulence in Couette- Taylor 
experiments have been discussed by Beck et al in p5f for four typical values of the Reynolds 
number Re and millions of experimental distances r (in units of the Kolmogorov length if). 
The data they present in their Fig. 3 can be organized as shown in the present Fig. 6, with 
q eq (oo) ~ 1.45, and the role played by the exponent 0.37 is that of q mix in Eq. (12). Inciden- 
tally, 0.37 is the value used by the Arimitsu's to produce (in one of their calculations) 
such distributions of velocity differences. It is worthy mentioning also that the value 1.45 is 
very close to 3/2 recently advanced by Beck for Lagrangian turbulence. 

Saddle-point dynamics of the Henon-Heiles system have been discussed by Soares et al 
33 1. In their Table 1 they show the numerical values obtained for 7 = l/{q eq — 1) as a 



function of the control parameter k, known to play a role analogous to a Reynolds number, 
hence k characterizes W and the k —>■ 00 limit corresponds to infinitely fine graining. The 
data presented in that Table can be organized as shown in the present Fig. 7, thus giving 
support to relation (12), with q eg (oo) ~ 2.81, and the role played by the exponent 0.35 is 
that of q mix . 

Summarizing, the basic picture which emerges from all the above considerations is as 
follows. If the system is strongly chaotic (in the sense that its spectrum of Lyapunov 
exponents includes positive values, i.e., exponential mixing), then the measure of ignorance 
(entropy) to be used is, as well known, Si, from which the entire Boltzmann-Gibbs statistical 
mechanics can be derived. If, however, the system is only weakly chaotic (zero Lyapunov 
exponent spectrum, and power-law mixing), then several indications exist which suggest 
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that we should instead use S qmix in what concerns a finite entropy production, and S qeq to 
discuss the corresponding equilibrium thermodynamics, both q mix and q eg being uniquely 
determined by the dynamics of the specific physical system (to be more precise, the graining 
degree, characterized here by W, also enters in the determination of q eq (W)\ it is the value 
q eq (oo) the one which only depends on the fundamental dynamics). As we see, in this scene 
there are two protagonists, namely mixing and equilibration. Although deeply entangled, 
they are different concepts, typically represented by two different values of q, one of them 
{Qmix) not above unity, and the other one {q e q) not below unity. These two values of q merge 
on a single one (q — 1) for Boltzmann-Gibbs thermostatistics (thus transforming power laws 
in exponentials). It is certainly allowed to think that this is perhaps at the origin of not few 
of the warm controversies in our community about the foundations of statistical mechanics! 

One of us (CT) acknowledges very fruitful discussions with E.G.D. Cohen, E.M.F. Cu- 
rado, A.C.N, de Magalhaes, A. Plastino, A.R. Plastino (see Appendix), D. Prato, I.D. Soares, 
H.L. Swinney, V.V. Uchaikin (see Appendix) and R.O. Vallejos. We are grateful to C. Beck 
et al for kindly providing to us the numerical data corresponding to Fig. 3 of |35|| . Finally, 
unforgettable hospitality from the main organizers of the meeting, namely P. Quarati, M. 
Lissia and A. Rapisarda, is also warmly acknowledged. This work has benefited from partial 
support from CNPq, PRONEX, CAPES and FAPERJ (Brazilian agencies). 

APPENDIX 

This Appendix focuses on an interesting point that was raised during the meeting by 
V.V. Uchaikin, and reports A.R. Plastino' s related remarks. 

Let us first remind, along the lines of |56| (see also f57fl), the q = 1 case. Let us assume 
that we have a continuous probability distribution p(x) (with J dx p(x) = 1, x being a one- 
dimensional real variable) and, using the entropic functional S*i[p(x)] — — j dx p(x) lnp(x), 
we wish to consider its discretization. In other words, we consider pi = p(xi)A, where 
A represents a graining for x, and i = 1,2, ...,W. It follows that, in the limit A — > 0, 
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5 1 [p(x)]~-Ell 1 p(^)ln £ f 



J2i=iP( x i) ]np( x i) + In A, hence, 



S 1 \p{x)]~S 1 ({p(x i )})+hiA. 



(13) 



We see that A — > leads to a divergence. This divergence is not surprising and corresponds 
|56| to the fact that it is necessary an infinite number of yes/no answers to resolve the 
uncertainty associated with a continuous distribution. In practice, it is chosen A = 1, 
therefore the continuous and discrete versions of the entropy provide the same result. The 
choice A = 1 corresponds, when we consider the passage from quantum to classical statistical 
mechanics, to the measure of phase space in units of h for every pair of conjugate mechanical 
variables. 

Let us now address this point for arbitrary q. Using the entropic functional S q [p(x)] = 



where we have used Eq. (1). Of course, this equation recovers Eq. (13) for q — 1. We 
see that, on top of an extra additive term, like in the q = 1 case, we have here an extra 
multiplicative term. Nevertheless, like in the q = 1 case, the choice A = 1 makes the 
continuous and discrete versions of the nonextensive entropy to coincide. 



l-J dx [p(x)]i 
9-1 



we straightforwardly obtain 



5 9 [p(x)]~A 1 -%({p(x i )})+ln 9 A, 



(14) 
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TABLE I. Historical steps of the foundations of BG statistical mechanics (both equilibrium 
distribution and entropic form), and their g / 1 counterparts. Gibbs 1902 and Kubo et al 1988 
refer not to the dates when the original works were essentially done, but rather to the books where 
they are reproduced. 
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FIG. 1. Traditional view (schematical) of the place of statistical mechanics for classical, quan- 
tum or relativistic dynamical systems with a finite or infinite number N of particles. Thermal equi- 
librium only occurs for conservative systems and is necessarily of the BG class, independently of the 
ordering of limits such as the t — > oo and the TV — > oo ones. This is certainly the case of short-range 
interactions which present no delicate singularity at say the origin. The micro-canonical ensemble 
can be seen as the particular case of the canonical one when the temperature diverges; the canonical 
one can in turn be seen as the particular case of the grand-canonical ensemble when all chemical 
potentials vanish. KSLPK stands for Kolmogorov-Sinai-Lyapunov-Pesin-Krylov, thus meaning the 
region where, for the particular case of classical dynamical systems, the Kolmogorov- Sinai entropy 
is positive, the Lyapunov spectrum includes a positive branch, the Pesin identity is nontrivial, 
i.e. connecting nonvanishing quantities, and Krylov's emphasis on exponential mixing being essen- 
tial in the foundations of BG statistical mechanics. The Sinai-Ruelle-Bowen (SRB) distributions 
are stationary states which can be both at or out from thermal equilibrium; the BG equilibrium 
distributions are particular cases. 
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PRESENT VIEW 



dynamical 
attractors 



arbitrary currents 





geometrizable 
dynamics 



^-region 

stationary states / closed systems 

thermal equilibrium (V q) 



FIG. 2. Present view (schematical) of the place of statistical mechanics for classical, quantum 
or relativistic dynamical systems with a finite or infinite number N of particles. The traditional 
ensembles are enlarged in the sense that q can differ from unity. Distributions similar to those 
occurring at thermal ^-equilibrium, or metaequilibrium (in the sense of metastability) , can occur 
even for dissipative systems. Ordering of limits such as the t — > oo and the N — > oo ones can 
be very relevant. For example, for conservative many-body systems including long-range interac- 
tions, the limt^oo limAr^oo ordering, which is the physical one, is non-Boltzmannian, whereas the 
limTv^oo hm^oo, physically unobservable, corresponds to BG statistics. The (/-region includes the 
1-region, which precisely is the KSLPK region of Fig. 1. 
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FIG. 3. Comparison between q- and stretched exponentials, (a) The circles have been calculated 
with a stretched exponential, and have been fitted with a g-exponential. (b) The circles have been 
calculated with a g-exponential, and have been fitted with a stretched exponential. The numerical 
discrepancies emerge only for x « 1 and for x » 1. 
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0.1 1 10 100 1000 

FIG. 4. Time evolution of the sensitivity to the initial conditions £. The crossover occurring 
in the limit Ai — > becomes apparent: the smaller Ai is compared with X qmix , the larger is the 
domain of validity of the power law £ oc The crossover time t* does not depend on Ai; 

the crossover time t** depends on Ai and diverges when Ai vanishes. 
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FIG. 5. Values of q eq obtained by Bediaga et al through the analysis of distributions of hadronic 
transverse momenta in electron-positron experiments. We have chosen the abscissa in such a way 
as to produce a linear form. R 2 is the square linear correlation factor. 



26 



1.20 

q 

1.16 



1.12 



1.08 



1.04 



1.00 




69 000 
133 000 
266 000 
540 000 



0.2 

x = {[ln(r/T])]/[(ln Re) 7/4 ]} 



0.4 0.6 

0.37 



FIG. 6. Values of q eq obtained by Beck et al through the analysis of distributions of velocity 
differences in Couette- Taylor experiments. We have scaled \n{r/rj) with (lni?e) 7 / 4 in order to 
produce data collapse for different values of Re. Also, we have chosen the exponent 0.37 in such a 
way as to produce a linear form for intermediate distances. The error bar on this exponent is of 
the order of 0.1 . 
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FIG. 7. Values of q eq obtained by Soares et al through the analysis of distributions emerging 
in the analysis of saddle point dynamics of the Henon-Heiles system. We have chosen the abscissa 
in such a way as to produce a linear form. R 2 is the square linear correlation factor. For the value 
A = 1/3 see [44]. 
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